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der conditions for electromagnetically induced transparency (EIT) . We find that the 
sharp edged axial density profile of an interacting condensate (due to a balance be- 
1/-^ I tween external trap and repulsive atomic interaction) is advantageous for obtaining 

CN ' ultra slow averaged group velocities. Our results are in good quantitative agreement 

. with a recent experiment report [Nature 397, 594 (1999)]. 
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1 Introduction 



Recent reports [1,2] of stopping a near resonant laser pulse inside highly ab- 
sorptive atomic vapors has attracted tremendous public interests. Physically, 
this magic like 'light trapping' is due to delicate quantum interference un- 
der the condition of electromagnetically induced transparency (EIT) [3-5]. In 
these experiments, a dynamical controlled pump field effects the conversion of 
a Raman-matched probe field into atomic Raman coherence of the medium. 
Broadly interpreted, it provides a mechanism for storage of quantum infor- 
mation encoded in the probe field. The experimenters also demonstrated the 
read-out by back-converting stored atomic Raman coherence into a new light 
field. EIT typically occurs under conditions of complete opaqueness which pre- 
vents the propagation of incoherent light; while a coherent propagation can 
succeed because of interference induced cancellation of one-path absorption. In 
several earlier experiments, ultra slow light propagation was observed in Bose- 
Einstein condensate (BEC) of Na atoms [~ 17 (m/s)[6] and ~ 1 (m/s)[7]], in 
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hot rubidium vapor [~ 90 m/s [8] and ~ 8 m/s [9]], and in a Pr:YSO crystal 
[~ 45 m/s [10]]. As the net reduction of hght speed is about the same order of 
magnitude in different host media, one may argue that being bose condensed 
is perhaps not so cruciaL Prom a practical point of view, however, the long 
coherence time of BEC does prove to be advantageous. Compared with a ho- 
mogeneous media [2,8-10], using a trapped atomic BEC as the EIT medium 
raises an important complication from its spatially inhomogeneous density 
profile (due to both external trapping and atom-atom interaction) [1,6,7]. A 
consistent interpretation therefore requires careful spatial averaging and theo- 
retical modeling. In this article we elucidate the effect of such a spatial average 
which also takes into account atom-atom interaction. 

Slow light propagation can be phenomenologically understood in terms of a 
high index of refraction (slow phase velocity) or a high dispersion-like feature 
(slow group velocity). It has been extensively studied for atomic media ex- 
hibit EIT [3]. Under typical conditions, such media are highly sought-after as 
there exist many potential technological applications, e.g. optical delay lines 
[6], quantum entanglement of slow photons [12], non-classical (e.g. squeezed) 
and entangled atomic ensembles [13], and quantum memories [1,2,5]. Other 
potential fundamental applications include: high nonlinear coupling between 
weak fields [6,11], quantum non-demolishing measurements and high precision 
spectroscopy [14], and as narrow-band sources for non-classical radiation [15]. 
In addition to a complete stoppage of light [5] , possibilities of even negative 
group velocities have been discussed [16]. We expect our investigation will 
shed light onto the realization of these applications in addition to providing a 
satisfactory understanding of slow light propagation in trapped atomic BEC. 

The paper is organized as follows: In Sec. II we briefly review our formulation of 
near-resonant light propagation in dispersive medium. Section III is devoted to 
the discussion of laser pulse propagation under conditions of EIT from a pair of 
Raman-matched pulses. In Sec. IV we discuss the widely used two-component 
model for the density profile of a trapped interacting BEC. In Sec. V we 
present numerical studies and compare them with the earlier experimental 
report of Hau et al. [6]. Finally we conclude in Sec. VI. 



2 Group velocity in a dispersive media 

Propagation of light in a non-magnetic, charge-free, dispersive medium is gov- 
erned by the second order wave equation 
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where E is the coarse grained electric field of the propagating hght and P is 
the induced macroscopic medium polarization. Limiting our discussion to a 
medium where spatial changes of the dielectric function are negligible within 
a wavelength of propagating light, the last term in the Ihs of Eq. (1) can be 
neglected. Although this is a standard approximation, we cautiously point 
out that it may lead to significant errors near the boundary of an interacting 
atomic BEC, where the gas density profile changes rapidly [17]. In our numer- 
ical simulation we also neglected this term based on an additional argument. 
For a transverse field propagating along z-axis, the coupling of V ■ P{r,t) is 
mainly along the orthogonal direction from the propagation direction. There- 
fore we only expect slight modifications to field distribution in the immediate 
neighborhood of the BEC boundary. We reconstruct the behavior of injected 
pulse in this region from an interpolation between its behaviors in and out 
of the condensate. The propagation is always started away from the edge of 
the BEC and the thermal component density contribution of trapped atomic 
gas in fully included. In the paraxial approximation limit, we make further 
simplification of Eq. (1) by neglecting the transverse Laplacian V'^E{f,t). 
This latter approximation is checked for self-consistency in comparison with 
experiments performed with a trapped BEC. 

— * 

We consider the propagation of a light pulse with central wave vector kg along 
the propagation (z) axis and carrier frequency ujq. In a dispersive medium, 
characterized by a complex refractive index n{k, u) at frequency uj and wave 
number k, the relation between the carrier frequency and wavevector is the 
dispersion relation n{kQ,u!Q)cuo — cko. Typically n(A;o,<^o) ~ 1 at resonance, 
therefore ouo — cko remains approximately valid throughout such a medium 
for near-resonant carrier frequencies. Employing the slowly varying phase 
and amplitude approximation for the complex amplitudes {V,S}, defined by 
{P, E} = ^[{V, £} exp {ikoz — iuJot)], we obtain 



dS Ids AndV ^ ., ^ 

The medium polarization is to be calculated from the linear response of atoms 
to the weak probe field E. It is customarily described in terms of an electric 
susceptibility x(r, t) according to 



V{z,t)^ j dz' j dt'x{z- z',t-t')S{z',t'). (3) 

All macroscopic physical quantities used are assumed to be derived from a 
coarse graining procedure where the medium is assumed isotropic microscopi- 
cally. According to the convolution theorem, the integral in the rhs of Eq. (3) 
can be expressed as. 
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V{z,t)^ j duj j dkx{k,uj)E{k,ujy^'-^\ (4) 

which is equivalent to V{k,uj) = x{k,uj)S{k,uj). To simphfy the notation, we 
will reference {k,uj) with respect to the origin (/cqjO^o)- Expanding x{k,uj) to 
first order around Xo = xl^o^'-^o) one obtains [18] 



which upon substituting into Eq. (2) yields [8,16] 
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Under the assumption of slow spatial variations we can neglect the spatial 
dispersion term ko{dx/ dk)Q. We now introduce a complex valued function 
Ng — 1 + 27rx(0) + 27ruJodx/duJo, whose real (imaginary) part will be called 
group index (phase correlation index). The real (imaginary) part of l + 27rx(0) 
will be called refractive (loss or gain) index. Unless otherwise stated, we use 
(.)' and (.)" to denote real and imaginary part of (.) throughout this article. 
The Wave equation (6) then simplifies to 



dz Vg dt 



(7) 



where the complex parameter a = 27riA;oX(0) mainly governs pulse attenuation 
or amplification while the complex velocity function Vg = c/Ng is mostly 
responsible for the propagation. In order to appreciate the kinematic meaning 
of Vg, it is helpful to cast the complex wave equation (7) into two coupled 
equations for real functions (p (phase) and U (amplitude) defined according to 



The propagation of a signal as well as the signal speed can only be defined for 
such real observables. When N'J pa 0, these two equations become uncoupled 
with Eq. (8) governing pulse propagation kinematics. Consequently Vg = c/N'g 
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Fig. 1. The group velocity for a pulse propagating through a three-level Raman EIT 
medium of a uniform density calculated according to Vg = ^{c/Ng). See later for the 
discussion gas density ~ 3.3x 10^^(cm-3), Aq = 589(nm), T31 = 0.67, Tsi = (27r)10^ 
(Hz), 9, = 0.567, and 7 = (27r) 10.01 (MHz). 




Fig. 2. Same as in Fig. 1 except we used Vg = c/^{Ng). 

can be defined as group velocity. Ambiguity arises when Eqs. (8) and (9) are 
coupled for N'^ 7^ 0. Our choice of group velocity definition Vg is equivalent 
to Vg = c/^{Ng), which is not the same as Vg = ^{du/dk) = K(c/iVg) [16]. 
When Ng is not negligible these two definitions generally lead to different 
results as can be exemplified by selected results presented in Figs. 1 and 2, 
which compares calculated group velocities according to the above two differ- 
ent definitions. The calculation of A^^ as well as the details of our model will 
be given in the next section. At zero probe field detuning N'J vanishes and 
the above two definitions become identical. For detunings where N'^ vanishes, 
however, the two results are drastically different as such points are singulari- 
ties for Vg = c/N'g. Thus, it seems more appropriate to treat the group velocity 
as a complex function in general, and assign a real physical meaning opera- 
tionally based on the pulse delays. The imaginary part causes non-trivial 
coupling between the amplitude and phase of a propagating pulse. As a re- 
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suit of this coupling, different propagation and transmission characteristics 
can occur. This situation is similar to the description of quantum tunneling in 
terms of a complex time and velocity [19]. When Ng is large, the amplitude 
and phase equations are strongly coupled, resulting in strong phase-amplitude 
correlations. This equation set can also be compared to the semi-classical laser 
equation [4] where similar situations arise. One important consequence of such 
amplitude-phase correlations is a 'false' gain phenomenon, or a perceived su- 
perluminal propagation [20]. To illustrate such effects explicitly, we consider 
a gas medium of uniform density. In this case, Eq. (7) can be solved with an 
ansatz exp [a{z — Zo)]f[t — {z — Zq)/v\, where Zq is the injection location of 
the incoming pulse and /(.) is an arbitrary complex valued function at the 
complex retarded time due to complex velocity v = C / {N'^ + iNg). Assuming a 
Gaussian temporal profile for the injected pulse at Zq as exp (— at^), the phase 
and amplitude functions U and can be solved 
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One can verify their correctness by direct substitution into Eqs. (8) and (9). 
We note that a non-vanishing A^^' contributes to U as an amplification. Usu- 
ally appreciable values of Ny are found in the anomalous dispersion regions 
where iV^ < 0. A negative A"^ results in a negative delay time (or advance) 
of pulse propagation, i.e. a superluminal propagation. Hence within our gen- 
eral model, superluminal propagation is accompanied by an amplification fac- 
tor exp [a{Ng{z — zo)/c)'^] in addition to the nominal loss/gain factor of the 
medium exp [a'{z — Zq)]. This suggests that even in an absorptive medium 
with a' — a superluminal pulse can tunnel through provided that the 

medium produces strong phase- amplitude correlation and is of a sufficient 
length. Essentially one needs a medium of length 



to beat the absorption factor. To illustrate this point, wc take a' = —0.1 (l//xm), 
A^ = —0.03c, A^' = 0.03c, typical numbers in the anomalous dispersion re- 
gion for near resonant pulse with a detuning A Ki O.27, and a — 0.44 (l//is)^, 
corresponding to a pulse of temporal width of 2.5 (/^s). We obtain a critical 
value oi L > Lc — 250 (//m) for superluminal propagation. 

In Fig. 4, we notice that A^ has a local minimum at A pa O.37. This local 
minima takes a negative value which signals a maximum negative delay time. 
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Fig. 3. Refractive index (solid curve) and loss index (dashed curve) for an atomic 
cloud with a density 3.3 x 10^^ (cm^"') corresponding the peak density of a trapped 
BEC at T ^ 450 (nK), slightly above Tc [6]. Other parameters used are: Aq = 589 
(nm), Tai = O.57, T21 = (27r)103 (Hz), Q. = O.567, 7 = (27r)10.01 (MHz). 
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Fig. 4. Group index (solid line) and phase correlation index (dashed line) (scaled 
by 1/c) under the same conditions as Fig. 3. 



However, in this case N'J is vanishingly small and it requires L — > cxd to beat 
the absorption. With careful analysis, one can choose an optimum value of 
Ng to take advantage of the amplification due to a nonzero Ng . A typical 
simulation for the physical observable U'^{z^t) is given in Fig. 5. For the pa- 
rameters used, a temporal advancement of the pulse, in other words a negative 
delay of 7.5 (/^.m) is found, which corresponds to a velocity higher than c by 
~ 30(m/s). It is possible to get higher super luminal velocities by decreasing 
the atomic vapor density. For instance, at p = 3.3 x 10^°(cm~^), which is 100 
times smaller than used in Figs. 3 and 5, the relevant parameters become 
iV^ = -3 X IQ-^c, A/"^' = 3 X IQ-^c, and a' = -1000(l/cm). In this case, super- 
luminal tunneling occurs for L > 2 (cm). These parameters are experimentally 
accessible and in fact are close to reported values in a recent experiment [21], 
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Fig. 5. Superluminal pulse propagation, or tunneling through an EIT medium. The 
negative delay time of 7.5(/xs) corresponds to a velocity higher than c by 30 (m/s). 

where N'g = —8 x 10~^c is measured in a Rb vapor of length L = 5 (cm). 

To develop an operational theory for the slow light propagation as in the exper- 
iment [6] , we need to include the above discussed amplitude-phase correlation 
as well as the non-uniform atomic gas density distribution. Furthermore, it 
would be necessary to include transverse diffraction losses as they might con- 
tribute significantly as the axial length of a trapped BEC increases. With all 
these complications, we have to resort to numerical simulations with available 
experimental parameters. Thus, we cannot simply define a group velocity for 
a pulse as c/Ng or some other form. We should define a viable speed in accor- 
dance with the actual propagation of a physical pulse. In this article we follow 
the operational definition as used in experiments [6,7]. We will concentrate 
on ultraslow pulse propagation rather than the equally likely superluminal 
effects. 

The existence of amplitude-phase correlation calls for caution in interpreting 
both slow and fast light propagation in a dispersive or active medium. Straight 
forward conclusions based on a linear response theory and simple definitions 
of group velocities can not be taken very seriously, especially when conflict- 
ing view with non-locality properties of electromagnetic fields arise. Strictly 
speaking, a pulse is composed of Fourier superposition of infinite wave trains 
with potentially different phase velocities. It is meaningful only when there 
exists a finite bandwidth around central frequency so that we can write [22] 



ko+5k 

E{z,t)= J E{k)e'^'-''^^^^'dk. (13) 

fco— <5fc 

Within a small bandwidth, linearization of the frequency around the center 
carrier gives uj{k) ~ ujq -\- Vg{k — ko) with Vg = dw/dk (evaluated at ko) the 
perceived group velocity. Such a velocity represents the speed of propagation of 
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constant pulse amplitude surface, a non-existent situation when strong phase- 
amplitude correlation exists. The large cross phase- amplitude due to loss/gain 
index (or the imaginary part of the group velocity) causes ambiguity since the 
constant amplitude surface follows a trajectory now in the complex plane. 
As a result, we recommend treating the group velocity as a complex valued 
function. Further discussions will depend on specifics of medium's electric 
susceptibility. This is the subject of the next section. 



3 A Raman- matched EIT medium of cold atoms 



Our study in this article is based on a medium of three-level A-type atoms. The 
excited level |3) is assumed to couple to the final state level |2) by a dressing 
laser beam with a Rabi frequency Q, and to the ground state |1) by the probe 
beam, whose propagation is the central issue of study. The dressing pump 
laser is assumed to be nearly Raman- matched with the probe laser (of carrier 
frequency ujq), which is detuned by A = ujq — uj:^ from the |3) |1) transition 
frequency cu^i. Initially all atoms are assumed to be in the ground state level 
In the weak probe field limit, a linear response calculation can be carried 
as only a small fraction of atoms are pumped out of their initial state. In this 
study we use the approximation that the initial density profile pi(r) being 
unchanged. Solving the density matrix equation for atoms, we determine the 
linear susceptibility to be x = rjXi, directly related to the steady state value 
for the |3) — > |1) matrix element. Here 77(f) = (3A|j^/327r^)pi(f) and A31 is the 
transition wavelength and [3] 



Xi(A,Q) = 




02 ^ ^ -1 



4r3i(r2i-2A)^ 



(14) 



Fsi is the decay rate from |3) to |1) and is the decoherence rate between 
1 2) and Typical detuning dependence of refractive and loss indices calcu- 
lated from this susceptibility are shown in Fig. 3. The twin peaks in the loss 
index are from dressed absorption fines (Autler-Townes doublet [23]), and the 
narrow valley between them is associated with the normal dispersion region in 
the refractive index curve. Within the valley region, absorptive loss is vanish- 
ingly small. Far away from this lossless normal dispersion region both curves 
diminish resulting in a marked increase of transmission, except around the 
absorption peaks, where optical response is dominated by normal dispersion. 
Anomalous dispersion regions are accompanied by significant absorptive loss. 
For comparison, the group and phase correlation indices are illustrated in Fig. 
4. In both normal and anomalous dispersion regions, we note that Ng and N'J 
are now dominated by the term (dx/dA) due to sharp variations in refractive 
and loss indices. In the normal dispersion region Ng > 0, and particularly in 
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the almost lossless region it attains a maximum on resonance,. This produces 
a minimal and well defined group velocity as N" ~ in the same location. In 
the anomalous region becomes negative or zero, corresponding to vanish- 
ingly small or negative delay times of propagation, as recently observed in [9]. 
The frequency derivative dx/dA is found to be rjdxi/dA with 



dA 



1 



Tai 4r3i(r2i-/A)2 



S2' 



4r3i(r2i-iA), 



(15) 



We can then define an auxiliary function h{A, Q) = Xi + (^ + <^3i)^Xi/^^) to 
express Ng = l + 27rrih{A, Q). We note that is maximum on resonance, and 
the trapped atomic density profile peaks at r = 0. Therefore an estimate of 
the minimum for group speed is Vg = c/[l + 27r?7(0)/i'(0, ^2)]. For a sample size 
of length L, we introduce dimensionless variables z — 2z/L and t — {Vg/L)t, 
and a normalized density profile function f{r)— pi(r)/pi(0). In terms of these 
parameters, the wave equation Eq. (2) becomes 



d£_ 
dz 
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< h{AM) 
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In a typical experiment with A = and F^i ^ , we find 



/n o\ •4r3ir2i 

'dxA £31 



(17) 



Thus, A^^' is indeed neghgible and we have Ng 27ra;3i77(4r3i/Q^). In this 
well defined group velocity exists and is given by 



Vg = C N' ^ -.^ CUa = ^7^. (18) 

For convenience, the last expression is in SI units {d^i/fi — 3A|ir3i/327r^). 
This result agrees with the one used for estimations of experimental results 
[6,25]. It is also the same as obtained earlier [11,26-28]. When applied to 
trapped atomic BEC at ultra low temperatures, inhomogeneous density profile 
needs to be included properly. When local density approximation is assumed, 
we expect the local light velocity should display rapid slowing down upon 
entering the BEC region from the thermal component and rapid acceleration 
when exiting it. This in general will result in different averaged velocity from 
simple estimations based on a uniform density. In order to elucidate this effect, 
it is necessary to specify an explicit form for the density profile of the BEC. 
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4 Spatial density profile of an interacting BEC 



In this section, we review briefly a simple analytical model for the density 
proflle of an interacting BEC as developed in Ref . [29] . We take the bose gas 
(below Tc) as composed of two components; a condensate part, whose density is 
computed using the Thomas-Fermi approximation and a background thermal 
atom part whose density is computed semi-classically assuming a continuum 
density of states in a harmonic trap [29,30]. The total ground state density is 
then found to be 



Pi« = ^9(/.-mrc-r) + *£5^|^, (19) 

where U = Anh'^ a^c / M . a^c is the atomic scattering length and 6{.) is the Heav- 
iside step function, gn{x) = J2j /f \ is the thermal de Broglie wavelength 
and P = 1/ {ksT). The external trapping potential is V{f) = Muj'^{r^ + e^z^)/2 
with cUr the radial trap frequency and e the aspect ratio. M is atomic mass. 
The chemical potential n is determined by the normalization N = J drpF{r) 
with N the total number of atoms. At high temperatures it is determined by 
solving Li3(z) = (T/Tc)^^C(3) for fugacity z = e^^. Li3(.) and (^(3) are the 
third order polylogarithm and Riemann-Zeta functions respectively. At low 
temperatures, it is found that [29] , 



N \ 2/5 



with htf the Thomas- Fermi approximation to chemical potential and the 
condensate fraction 



N 



1 - 



T_\^ _ C(2) / T 
Tc) ^C(3) \Tc 



2 r 



1 - 



^ \ 3' 

7b 



1 2/5 



(21) 



with a scaling parameter s. 



^ = \m'Hl^N^"^) . (22) 



Oho = V ^/ {Me^/^ijJr) denotes the average harmonic oscillator length scale [31]. 
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Temoerature fnK) 

Fig. 6. Comparison of an interacting (solid line) and ideal (dashed line) models for 
the temperature dependence of operationally defined group velocity. We have used 
the same parameters as in in Fig. 3. In addition, AT = 8.3 x 10^, Qg = 2.75 (nm), 
LUr = (27r)69 (Hz), co^ = (27r)21(Hz), M = 23 (amu) as for the experiment [6]. 
i^o = 15 (/xm). 

5 Numerical results and discussion 

In this section, we discuss several aspects of our numerical investigation for 
the experiment [6]. 

5.1 Group velocity in ideal and interacting BEC 

According to experimental prescriptions [6,7] we take the operational ap- 
proach, first determine the delay time according to 




where R is the pinhole radius introduced to selectively detect propagated light 
from the near axial region of the atomic cloud. As justified before [24] we will 
ignore the spatial {R) dependence of L due to the shape of the cloud for small 
pinhole sizes and relatively long optical path length L [25]. Within the local 
density approximation, the local group velocity is given by Vg{f) = c/N'g, and 
its average is then calculated according to {vg) = L{T)/td{T) with L{T) is 
the on-axis cloud size at a given temperature T. In Fig. 6 we compare the 
calculated results for an interacting and an ideal atomic gas. For the ideal gas 
case, our results simply reproduced those in Ref. [25]. Whenever possible we 
have used parameter values taken from Hau's experiment [6]. At high tem- 
peratures both models are in agreement with the experimental data. At low 
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T(nK) 



Fig. 7. Group velocity is calculated with 

interacting BEG for Osc = 7, 5.75, 3.75, l(nm) corresponding to solid, dot-dashed, 
dashed, and dotted curves, respectively. Other parameters are the same with those 
in Fig. 6. 

temperatures, our model including atom-atom interaction gives results of the 
correct order of magnitude, while those based on ideal non-interacting atoms 
[25] are two orders of magnitude too low. We also note that the sharp low 
temperature dependence near the BEC transition temperature Tc becomes 
significantly smoother, i.e. closer to the experimental observations when the 
atom interaction is included. Subsequently the calculated (vg) does not fall 
down as rapidly at lower temperatures. Our model including atom-atom in- 
teraction also predicts that (vg) increases with the atomic scattering length 
as illustrate in Fig. 7, where the group velocity is calculated for vari- 
ous scattering length parameters Ogc = 7,5.75,3.75, 1 (nm), corresponding to 
s = 0.3982,0.368,0.31,0.183, respectively. We note that the two component 
density profile works particularly well for s ^ 1, but becomes less accurate 
when s > 0.3 [29]. Thus, we conclude that the inclusion of interatomic inter- 
actions improves upon the ideal gas model prediction both qualitatively and 
quantitatively. In fact, our results with atom interaction give the correct or- 
der of magnitude as reported in the experiment [6] . Although we are still not 
near the reported value v ~ 17 (m/s) for T ~ 50 (nK). Inclusion of condensate 
number fluctuations, evaporation dynamics, and a self-consistent treatment of 
density profile might bring in further improvements. 



5.2 Axial propagation of probe pulse 



At temperatures well below Tc, both axial and radial density profiles of an 
atomic cloud are dominated by its condensate component. For the harmonic 
trapping potentials considered here, the Thomas-Fermi approximation results 
in an inverted paraboloid density profile. Consequently, the spatial behavior 
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Fig. 8. Axial profile of the scaled group index N'g/c at T = 43 (nK) and A = 0. 
Ng/c is very small with a peak value ~ 10~^^(s/m). All parameters are the same 
as in Fig. 6. 




Fig. 9. Propagation of a resonant pulse with A = 0, T = 43 (nK). All parameters are 
the same as in Fig. 6. Position and time are given in dimensionless units explained 
in the text. 

of the susceptibility in the linear response regime also follow such a behavior 
as shown in Fig. 8. A typical simulation for the time delay of a resonant pulse 
exp (— lOOt^) after its passage through such a medium is is shown in Fig. 9. 

In this simple wave packet propagation, we only considered axial propagation 
along z. A discussion of paraxial effects and diffraction losses are given later 
in the next subsection. This one dimensional approach is shown to be valid 
in several recent slow pulse experiments where the incident Gaussian beam 
radius is much larger than the pinhole radius and the effects of the far radial 
tail become negligible. We find that upon entering the medium, the probe 
pulse rapidly slows down due to the sharply increased refractive index of the 
BEG component. It then propagates with a minimum speed around the center 
of the cloud, finally accelerates back to its initial vacuum speed as it exits near 
the edge of the cloud. According to this physical picture, the average delay 
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position 



time 



Fig. 10. Absorption of an off-resonant pulse with detuning A = 87, 7 = (27r) 10.01 
(MHz), and T = 43 (nK). All other parameters are the same as in Fig. 6. Position 
and time are given in dimensionless units explained in the text. 

which is used in determining the average group velocity is then lower than the 
optimal minimum group velocity at the trap center. Indeed, from the figure 
we see that the piilse arrives at the exit edge at about 0.6 of the predicted 
delay time estimated from a uniform density. In Ref. [25], a group velocity of 
9 (m/s) was estimated for an interacting gas of effective uniform density with 
N = 10^ atoms, using the peak density from Thomas-Fermi approximation. 
When taking into account the condensate edge effects by directly propagating 
the wave equation in the scaled form, we obtained an improved estimation 
of 15 (m/s), which agrees quite well with the reported value of 17 (m/s). 
The average group velocity can be made very close to the desired minimum 
speed by choosing a density profile with sharp edges so that propagating light 
traverses the medium with the minimum speed for a longer fraction of total 
time spent in the medium. Ideally, a step function type density distribution 
along the axial direction would achieve this. 

To complete this section, we also briefly studied absorption of the pulse. We ex- 
amined the absorption length in detail using the propagation equation Eq.(16). 
At a large detuning of A = 87 we observe a penetration into only the first 15% 
of the cloud as demonstrated in Fig. 10. By controlling the absorption length 
via detuning, interaction of light with only a selected part of the condensate 
might be realized in such a simple way. 

5.3 Transverse diffraction effects 

The discussion of radial density profile effects requires the inclusion of prop- 

agational diffraction. In practice, one can always choose appropriate trap pa- 
rameters and suitably focused incident lasers to avoid most diffractive losses 
[6]. To examine the diffraction of a probe beam as it propagates through the 
condensate, we now include the previously neglected term to the wave 
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Fig. 11. Diffraction pattern of a beam of radius 0.5 (mm) as it propagates through 
a Bose condensate at T = O.STc- All other parameters are the same as in Fig. 6. 

equation. Employing a standard pseudo-spectral split operator method for 
wave packet propagation [32] we have simulated the propagation of a pulse 
of radial cross section size 0.5 (mm) in accordance with the experiment of 
Hau et al. [6]. Figure 11 illustrates its axial evolution of the time averaged 
intensity pattern. In general, we find essentially one dimensional propagation 
equation Eq. (16) adequately describes near axis propagation as performed in 
the experiment [6,7]. In the Fig. 11, we only plotted the y = plane because 
of its axial symmetry. Within the pinhole dimension set at ~ IS/^m by exper- 
imental arrangement, and for detection close to the atomic cloud, diffraction 
effects are seen to be qualitatively neghgible. Thus the large incident wave 
radius effectively assures a plane wave signal for the cloud size which is much 
smaller. The only significant effect of the transverse diffraction is then the 
diffraction loss of the pulse amplitude. In the experiment by Hau et al. [6], it 
was estimated that 80% transmission of the probe pulse should occur under 
the EIT condition, yet even at high temperatures the observed transmission 
was limited to about 30% [6]. Since Fig. 11 indicates a lower output intensity 
than found from the one dimensional axial propagation model, we suspect that 
the additional reduction in the pulse amplitude can be due to the transverse 
diffraction loss. 



6 Conclusion 

We have performed a detailed investigation of Raman-matched propagation 
of a probe laser pulse through a BEC. We have discussed of the interpreta- 
tions and definitions of the group velocity, and argued that a complex valued 
group velocity should be used when strong amplitude-phase correlation ex- 
ists. We have adopted the highly successful two component model for atomic 
condensate density profile and compared the differences of slow light propa- 
gation through an interacting and an ideal gas BEC. Based on the broadened 
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condensate profile of an repulsively interacting Na atom BEC, we have ob- 
tained group velocity values of the same order of magnitude as observed in 
the experiment [6]. Furthermore our simulations for this experiment explains 
the qualitative behavior observed in the temperature dependence of the group 
velocity at low temperatures [6] . In the ultracold regime deep below the BEC 
transition temperature, including the atom-atom interaction and the inhomo- 
geneous density profile results in two orders of magnitude improvement in the 
quantitative description of the group velocity behavior [25]. We also observed 
an increase in the group velocity at larger atomic scattering length. 

When comparing the detailed spatial averaging of operationally defined delay 
time, we find that a sharp edged density profile yields lower average group 
velocity (vg), i.e. closer to its local minimum Vg{0), typically obtained based on 
the estimation of group velocity using the peak density. Our detailed numerical 
simulations confirm (vg) ~ Vg{0)/0.6, which gives essentially the same (vg) 
value as reported [6]. 

We have also compared the One dimensional numerical simulations with a 
paraxial approximation including the diffraction effect during the probe pulse 
propagation. We find that under the experimental conditions, the use of a 
small pinhole to select probe light rays close to the x-axis is a reasonable 
approach as the diffractive effects are negligible apart from a total loss of the 
overall signal level. Thus, our studies also justify the one dimensional approach 
normally used in EIT studies. 



7 acknowledgments 

During the course of this work, one of us (L.Y.) was a participant of the recent 
workshop on quantum degenerate gases at the Lorentz Center, University of 
Leiden. L.Y. thanks H. Stoof for the hospitality. This work is supported by 
the NSF grant No. PHY-9722410. 



References 

[1] C. Liu, Z. Button, C. H. Behroozi, and L. Vestergaard Hau, Nature 409, 490 
(2001). 

[2] D. F. Phillips et al, Phys. Rev. Lett 86, 783 (2001). 

[3] S. E. Harris, Physics Today 50, 36 (1997); E. Arimondo, Progress in Optics 
(ed. E. Wolf), 257, (Elsevier Science, Amsterdam, 1996). 



17 



[4] M. O. Scully and M. S. Zubairy, Quantum Optics, (Cambridge University Press, 
1997). 

[5] M. Fleischhauer and M. D. Lukin, Phys. Rev. Lett. 84, 5094 (2000). 

[6] L. V. Hau, S. E. Harris, Z. Dutton, and C. H. Beliroozi, Nature 397, 594 (1999). 

[7] S. Inouye, R. F. Low, S. Gupta, T. Pfau, A. Gorlitz, T. L. Gustavson, D. E. 
Pritchard, and W. Ketterle, ( |cond-inat /0006455D . 

[8] M. M. Kash, V. A. Sautenkov, A. S. Zibrov, L. Hollberg, G. R. Welch, M. D. 
Lukin, Y. Rostovtsev, E. S. Pry, and M. O. Scully, Phys. Rev. Lett. 82, 5229 
(1999). 

[9] D. Budker, D. F. Kimball, S. M. Rochester, and V. V. Yashcuk, Phys. Rev. 
Lett. 83, 1767 (1999). 

[10] A. V. Turukhin, J. A. Musser, V. S. Sudarshanam, M. S. Shahriar, and P.R. 
Hemmer, ( quant-ph/0010009| ). 



[11 
[12 

[13; 

[14 
[15 

[16; 
[ir 

[is; 

[19 

[2o; 

[21 
[22 
[23; 
[24 
[25; 



S. E. Harris, J. E. Field, and A. Kasapi, Phys. Rev. A 46, R29 (1992). 

M. D. Lukin and A. Imamoglu, Phys. Rev. Lett. 84, 1419 (2000). 

M. D. Lukin, S. F. Yelin, and M. Fleischhauer, Phys. Rev. Lett. 84, 4232, 
(2000). 

A. B. Matsko, Y. V. Rostovtsev, H. Z. Cummins, and M. O. Scully, Phys. Rev. 
Lett. 84, 5752 (2000). 

M. Fleischhauer, M. D. Lukin, A. B. Matsko, and M. O. Scully, Phys. Rev. 
Lett. 84, 3558 (2000). 

O. Kocharovskaya, Y. Rostovtsev, and M. O. Scully, ( |quant-ph / 00010*5^ ) . 

F. Dalfovo, L. Pitaevskii, and S. Stringari, Phys. Rev. A 54, 4213 (1996). 

S. E. Harris, J. E. Field, and A. Kasapi, Phys. Rev. A 46, R29 (1992). 

R. Landauer and Th. Martin, Rev. Mod. Phys. 66, 217 (1994), and references 
therein. 

A. M. Steinberg, P. G. Kwiat, and R. Y. Chiao, Phys. Rev. Lett. 71, 708 (1993); 
A. Enders and G. Nimtz, Phys. Rev. E 48, 632 (1993); R. Y. Chiao and A. M. 
Steinberg, Prog. Optics XXXVII, 345 (1997). 

A. M. Akulshin, S. Barreiro, and A. Lezama, Phys. Rev. Lett. 83, 4277 (1999). 
J. A. Stratton, Electromagnetic Theory, (McGraw-Hill, N.Y., 1941). 
S. H. Autler and C. H. Townes, Phys. Rev. 100, 703 (1955). 
Ozgiir Miistecaplioglu and L. You, Phys. Rev. A (submitted, 2000). 

G. Morigi and G. Agarwal, Phys. Rev. A 62, 013801 (2000). 



18 



[26] R. Grobe, F. T. Hioe, and J. H. Eberly, Phys. Rev. Lett. 73, 3183 (1994). 

[27] M. Xiao, Y.-Q. Li, S.-Z. Jin, and J. Gea-Banacloche, Phys. Rev. Lett 74 666 
(1995). 

[28] A. Kasapi, M. Jain, G. Y. Yin, and S. E. Harris, Phys. Rev. Lett. 74 2447 
(1995). 

[29] M. Naraschewski and D. M. Stamper-Kurn, Phys. Rev. A 58, 2423 (1998). 

[30] V. Bagnato, D. E. Pritchard, and D. Klcppncr, Phys. Rev. A 35, 4354 (1987). 

[31] S. Giorgini, L. P. Pitaevskii, and S. Stringari, Jour. Low Temperature Phys., 
109, 309 (1997). 

[32] B. M. Garraway and K-A. Suominen, Rep. Prog. Phys. 58, 365 (19995). 



19 



